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We present deep near-infrared (J, H, Ks) images toward an embedded cluster 

which lies in a C^^O clump in the cloud C of the Vela Molecular Ridge. This 

^ . cluster has at least ~ 350 members and a radius of ~ 0.5 pc. The stellar surface 

^ . number density is approximately 3000 stars pc~^ in the central 0.1 pc x 0.1 

\Q . pc region of the cluster. This is much higher than most of the young clusters 

O. 

■^ . near-infrared excess fraction with those of other embedded clusters, we estimate 

that the age of this cluster is approximately 2-3 Myr. This cluster exhibits an 

p . excess of brighter stars in its central region, from which we conclude that the 

Q ! more massive stars are located near the cluster center. 

c^ . Subject headings: ISM: clouds — ISM: individual (Vela Molecular Ridge) 



stars: formation — stars: luminosity function, mass function 



1. INTRODUCTION 



The progress of NIR array detectors provided opportunities to observe stellar contents 
in molecular clouds. Many studies of these regions revealed the presence of stellar clusters 
which are deeply embedded in molecular clouds (Lada et al. 1991; Carpenter et al. 1993; 
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Hodapp 1994). Statistical studies revealed that the most of young stellar objects in giant 
molecular clouds (GMCs) are formed in embedded clusters (Lada & Lada 2003). These em- 
bedded clusters are typically found to have their most massive stars in the central region (the 
position whose stellar densities are the highest) (Zinnecker, McCaughrean, & Wilking 1993; 
Hillenbrand & Hartmann 1998). Several embedded clusters also shows mass segregation; 
more massive stars tend to be located in the central region [e.g., Trapezium (Hillenbrand 
& Hartmann 1998); NGC 2024, NGC 2071 (Lada et al. 1991); NGC 3603 (Niirnberger & 
Petr-Gotzens 2002)]. However, Carpenter et al. (1997) pointed out that there is no evidence 
of mass segregation in the Mon R2 cluster. Unfortunately, the mass segregation has been 
investigated only for a few embedded clusters. We need more samples for statistical studies. 

In this paper, we study an embedded cluster which lies in a C^^O clump identified by 
Yamaguchi et al. (1999) in the cloud C of the Vela Molecular Ridge (VMR). The VMR 
is a GMC complex, being composed of at least four GMCs (named A, B, C, and D) with 
individual masses exceeding 10^ M© (Murphy & May 1991). The clouds A, C, and D are 
located at 700 ± 200 pc, and the cloud B is located at ~ 2 kpc (Liseu et al. 1992). The 
cloud C is the richest in molecular gas and the least evolved of the complex (Yamaguchi et 
al. 1999). This embedded cluster was recently discovered by Massi, Lorenzetti, & Giannini 
(2003) in their NIR imaging survey toward Class I sources in VMR (Liseu et al. 1992; 
Lorenzetti, Spinpglo, & Liseu 1993). Massi, et al. (2003) found that high mass stars which 
excite the H II region RCW 36 lie in this cluster (hereafter RCW 36 cluster). However, due 
to the relatively small image dimensions, their observations were not able to cover the whole 
cluster. Therefore, the fundamental properties of the RCW 36 cluster, such as the number 
of the stars, radius and the stellar distribution, have been still unknown. Our observation 
first provides the detailed census of the RCW 36 cluster. 



2. OBSERVATIONS AND DATA REDUCTION 

The NIR images of RCW 36 were obtained on 2002 March 9 with the NIR camera SIR- 
lUS (Simultaneous three-color InfraRed Imager for Unbiased Survey) mounted on the IRSF 
(InfraRed Survey Facility) 1.4 m telescope of Nagoya University at Sutherland, South African 
Astronomical Observatory. SIRIUS is equipped with three 1024 x 1024 pixel HgCdTe arrays. 
Dichroic mirrors enable simultaneous observations in the J, if, and Ks bands (Nagashima et 
al. 1999; Nagayama et al. 2002). The image scale of the array is 0".45 pixel"^, giving a field 
of view of 7'. 7 x 7'. 7. Ten dithered frames were observed as a set of exposures. We observed 
three sets of the object and sky frames alternatively with a 30 second exposure. We also 
observed one set of the object and sky frames with a 5 second exposure for bright sources. 
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Seeing condition was ~ I'M (FWHM) in the Kg band. We observed the standard star 
9136 in the faint near-infrared standard star catalog of Persson et al. (1998) for photometric 
cahbration. 

We apphed the standard procedures of near-infrared array image reduction, including 
dark current subtraction, sky subtraction, and flat-fielding, using the IRAF (Imaging Re- 
duction and Analysis Facility) ^ software package. Identification and photometry of point 
sources were performed by using the DAOPHOT packages in IRAF. In order to avoid the 
affection for the photometry from the nearby stars, we used a radius of ~ 3 pixels (1".35) in 
the PSF fitting in the DAOPHOT task. 

The limiting magnitudes (at 10 a) are estimated to be ~ 19.0, 18.0, and 16.6 in the J, 
H, and Ks bands, respectively. A total number of 823 sources are detected in the Ks band 
with a photometric error of < 0.1 mag, and 610 sources are identified in the three bands with 
a photometric error of < 0.1 mag and with a positional agreement of < 0.5". In order to 
determine the completeness limit of the photometry, we added artificial stars to the central 
region of the image (dotted square in Fig. 2) with 17" spacing. We define the magnitudes at 
which 90% of stars were detected as the completeness limits. They are 18.2, 16.9, and 16.3 
at J, H, and Ks, respectively. 



3. RESULTS 

3.1. NIR image 

In Fig. 1, we show a composite JHKs color image of the RCW 36 cluster. It clearly 
shows that the stars concentrate at the center of the image. The brightest star located at 
8^59^27^5, -43°45'27" (J2000) is the exciting star (08 or two 09 stars; Verma et al. 1994) 
of RCW 36 (Massi, et al. 2003). At optical wavelengths, the stellar cluster including the 
exciting star is hidden by a dark-lane across RCW 36. The radio continuum peak position 
is ~ 35" west of the exciting star (shown by plus in Fig. 1; Walsh et al. 1998); this coincide 
spatially with a bright nebula. An H2O maser is detected ~ 30" west of the exciting star 
(shown by asterisk in Fig. 1; Braz & Epchtein 1982); there is no NIR counterpart. The 
existence of the H2O maser means that star formation is still ongoing in the RCW 36 cluster. 
There are several bright-rimmed clouds at ~ 1' west, ~ 1' south, and ~ 2' south-east of the 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science. 



-4- 



exciting star. They are likely to be dense cores which were revealed by evaporation of a 
molecular cloud due to UV photons from the exciting star. 



3.2. Stellar Census in the RCW 36 Cluster 

3.2.1. Stellar Surface Number Density 

In Fig. 2, we show the contour map of the stellar surface number density (SSND) derived 
from 823 sources detected at Ks with a photometric error of < 0.1 mag. The SSND is 
obtained by counting stars in a 30" x 30" area every 15" over the Ks image. At the distance 
of 700 pc, the unit cell of 30" x 30" corresponds to ~ 0.1 pc x 0.1 pc. The frequency 
distribution of the star counts (in stars per unit cell) can be fitted by the Poisson function 
with a mean of 2.0 stars per unit cell, plus the wing where excess of counts due to the central 
clustering appear. The SSND shows two peaks. The northern peak coincides spatially with 
the exciting star and shows the highest SSND. The southern one coincides spatially with one 
of the bright rims at ~ 1' south of the exciting star. 



3.2.2. Radial Profile 

In Fig. 3, we show the radial profile of the SSND for the RCW 36 cluster. The SSND 
is measured in circular annuli, centering on the location of the peak position of the SSND. 
Profiles are calculated using annuli with equal radial steps (=15"). The surface density 
profiles of embedded clusters can be fitted with 1/r profiles (e.g., McCaughrean & Stauffer 
1994; Jiang et al. 2003; Muench et al. 2003). They also can be fitted with King profiles (e.g., 
Horner, Lada, & Lada 1997; Niirnberger & Petr-Gotzens 2002; Teixeira et al. 2004). If we 
can neglect the tidal radius. King profiles can be written as 

/./ X /o 

fir) 



I + {r/rcY 

where /o is the central surface density and r^ is the core radius (King 1962). Ideally, the 
radial profile should become zero at r ^ oo. However, in the observation of the actual 
clusters, the radial profile does not become zero but a constant at r -^ oo due to the field 
stars which contribute as a background offset. Therefore, in order to take the contribution 
from the field stars into account, we added a constant term as a fitting parameter to fitting 
functions. 

The data of the RCW 36 cluster can be well fitted with King profiles better than with 
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1/r profiles; a departure from King profiles arising around r ~ 0.2 pc is due to the southern 
SSND peak. King profile fitting yields the central density of /o = 3200 stars pc~^ and the 
core radius of r^ = 0.08 pc. Since the column density of isothermal gas sphere in hydrostatic 
equilibrium has a radial dependence of r~^ (Yun & Clements 1991), the 1/r dependence of 
the SSND radial profile in the embedded clusters is likely to be a footprint of the density 
profile of their parental cloud core. On the other hand, King models describe solutions for 
stellar systems in dynamical equilibrium, such as globular clusters. The secondary peak of 
the SSND, which causes the departure around r ~ 0.2 pc in Fig. 3, might be a relic of the 
primordial substructure that has not fully merged into main body of the RCW 36 cluster 
(Scally & Clark 2002). Since King profiles better fits the observed radial profile than 1/r 
profiles overall except for the secondary peak, the dynamical equilibrium would have been 
almost, if not completely, established in the RCW 36 cluster. 



3.2.3. Determination of the Cluster Properties 

In order to obtain the fundamental properties of the RCW 36 cluster, such as the number 
of the stars and radius, we have to estimate the degree of the field star contamination. The 
frequency distribution of the star counts per unit cell (used to derive SSND in §3.2.1) can 
be fitted by the Poisson function with a mean of 2 stars per unit cell. This means that 
uniform background stars are distributed in the whole observed field. Hereafter, we take this 
value (2.0 stars/unit cell ~ 190 stars pc~^) as the background. This estimation is reasonable 
because the background value yielded by King profile fitting in the previous section is in 
good agreement with it. 

We can estimate the number of the RCW 36 cluster members by subtracting background 
from 823 sources detected at Kg with a photometric error of < 0.1 mag. The total expected 
background population in the observed field is 474 [= 2.0 x 4 (stars arcmin"^) x 7.7 x 7.7 
(area of observed field)]. Then, the total number of the cluster members is estimated to be 
349; we note that this estimation is a lower limit down to our limiting magnitude. The radial 
profile merges with the extended background at r ~ 0.5 pc (see Fig. 3). We take this as the 
cluster radius. Considering the background contamination, the SSND of the central 0.1 pc x 
0.1 pc region of the cluster is approximately 3000 stars pc~^; this is in good agreement with 
the central density yielded by King profile fitting in the previous section. This is different 
from the value of ~ 4800 stars pc~^ derived by Massi, et al. (2003), but we confirmed in our 
data that this inconsistency is due to their unit cell size (20" x 20"). In table 1, we show 
the comparison of the central SSND with other embedded clusters within 1 kpc of the Sun. 
The central SSDN of the RCW 36 cluster is much higher than most of the young clusters 



-6- 



but lower than the Trapezium cluster. 

The theory borrowed from the traditional theory of low-mass stars that stars are formed 
through gravitational collapse of dense molecular cores can only account for those stars with 
masses less than 10 M0, because radiation pressure of massive stars can halt the collapse 
and reverse the infall. Bonnell, Matthew, & Zinnecker (1998) presented result of numerical 
experiments that the massive (M > IOMq) stars can be formed through coalescence of two 
or more lower mass stars in the center of rich, dense stellar clusters. The most massive star 
in the RCW 36 cluster (08-09) is located at the cluster center; this might imply that it was 
formed through this mechanism. Hillenbrand (1995) suggested that there is a correlation 
of the cluster density with the maximum stellar mass in the cluster (i.e. more centrally 
condensed clusters have more massive stars). In fact, the most massive stars in R CrA, 
IC 348, and Mon R2 are less massive than that of the RCW 36 cluster, and Trapezium 
has a more massive star (see table 1). However, the most massive star of NGC 2024 is 
approximately the same as that of the RCW 36 cluster. We note that the central SSND of 
NGC 2024 might have been estimated to be too low, because the observation toward NGC 
2024 (Lada et al. 1991) was made with relatively small limiting magnitude (K ~ 14 mag) 
and relatively larger pixel scale (~ 1".3 pixel ~^). If we assume a distance of 415 pc and 
mean extinction of Ay ~ 10.5 mag toward NGC 2024 (Haisch, Lada, & Lada 2000), the 
limiting magnitude corresponds to an extinction-corrected absolute magnitude of M/^-^ ^^ ~ 
4.9 mag. For an assumed age of 1 Myr, their observation could detect down to ~ 0.09 M© 
(Baraffe et al. 1998). In order to take low-mass members of NGC 2024 into account down to 
our observation limit [~0.04 Mq assuming a distance of 700 pc, mean extinction of Ay ~ 8.1 
mag, and an age of ~ 3 Myr (see below)], we use the slope of Trapezium IMF (Muench et 
al. 2002). The estimated central SSND of NGC 2024 (~ 1700 stars pc'^) is still lower than 
that of RCW 36. However, the larger pixel scale might affect the estimation of the central 
SSND more seriously, because the radial profile of the SSND for NGC 2024 decreases in the 
cluster center (see Fig. 8 of Lada et al. 1991). 



3.3. color-color diagram 

In Fig. 4, we show the J-H v.s. H-K color-color diagram for the 610 sources identified 
in three bands (see §2). Their colors are transformed to the CIT system with the color 
equations^ (Y. Nakajima, et al. in preparation). The sources falling into a region of infrared 
excess (the right of the reddening band) in the color-color diagram are considered to be 



^http://www.z. phys.nagoya-u.ac.jp/~sirius/about/color _e.html 
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likely cluster members with optically thick, circumstellar disks. The total of 77/349 (22%) 
of the cluster members exhibit infrared excesses greater than their photometric errors (we 
note this excess fraction might be just a lower limit because several cluster members were 
not detected in three bands). We also estimate the mean extinction toward the RCW 36 
cluster to be Ay ~ 8.1 mag by de-reddening all sources in Fig. 4 back to this CTTS locus 
(Meyer, Calvet, & Hillenbrand 1997) along the reddening vector. 



3.4. K.S Luminosity Function 

In Fig. 5, we show the Ks band luminosity function (hereafter KLF) for the RCW 36 
cluster. In order to obtain the KLF for the cluster members (cluster KLF) , we must remove 
background contamination. We use the sources which are included in the regions whose 
SSND is smaller than the background plus la (off-cluster region) to construct the off-cluster 
KLF. The off-cluster region is approximately the same as the region r > 0.5 pc from the 
cluster center. Thus, most of the sources contained in the off-cluster region are supposed to 
be background stars (see Fig. 3). The off-cluster KLF is shown with open squares and dashed 
line in Fig. 5a. By subtracting the (scaled) off-cluster KLF from the KLF which includes 
both the cluster members and the field stars (raw KLF; shown in Fig. 5a with filled circle 
and solid line), we can obtain the cluster KLF (Fig. 5b). The cluster KLF has a turnover 
with a peak at Ks ~ 14 mag, and then decreases down to the completeness limit. 



4. DISCUSSION 

4.1. Age of the RCW 36 cluster 

NIR imaging surveys revealed that, in embedded clusters, the fraction of the sources 
which possess NIR excess decreases with the cluster age (Lada 1999; Haisch, Lada, & Lada 
2001b). The NIR fraction derived from JHK imaging is approximately 50 % for the young 
(< 1 Myr) embedded clusters [e.g., 50 ± 7% in Trapezium (Lada et al. 2000); 58 ± 7% in 
NGC 2024 (Haisch, Lada, & Lada 2000)]. The fraction decreases to approximately 20 % for 
more evolved (~ 2-3 Myr) embedded clusters [e.g., 21 ± 5 % in IC 348 (Haisch, Lada, & 
Lada 2001a); 16 ± 3 % in NGC 2316 (Teixeira et al. 2004)]. The NIR excess fraction of the 
RCW 36 cluster (22 %; see §3.3) implies that the age of the RCW 36 cluster is approximately 
2-3 Myr. 

A similar age is obtained by comparing the KLFs. The comparison of the KLFs of 
RCW 36, Trapezium and IC 348 is shown in Fig. 6. The KLFs for Trapezium and IC 348, 



by Muench et al. (2003), are corrected to absolute magnitudes. Reddening corrections were 
performed for each KLF by shifting to brighter magnitude by each mean extinction [Av ~ 
8.1 mag in RCW 36 (see §3.3), Av ~ 2.4 mag in Trapezium (Prosser et ah 1994), Av ~ 
4.5 mag in IC 348 (Lada & Lada 1995)]. Trapezium KLF and RCW 36 KLF are scaled to 
contain the same number of stars as IC 348. The shape of the KLF of RCW 36 is similar 
to that of Trapezium, but is shifted toward the fainter magnitude. On the other hand, the 
KLFs of IC 348 and RCW 36 show good agreement in shape and they have the same broad 
peak at similar absolute magnitudes. 

Muench et al. (2000) presented the result of numerical experiments that the KLFs for 
young clusters which have the same underlying IMF evolve in a systematic manner with 
increasing age; they evolve to fainter magnitudes. The IC 348 IMF is nearly identical to 
the Trapezium IMF (Muench et al. 2003), and they are similar to that derived for other 
embedded clusters, open clusters, and field stars (e.g., Lada & Lada 2003, and the references 
therein). Hence, if we assume a universal IMF for the RCW 36 cluster, we infer that the 
RCW 36 cluster has approximately the same age as IC348, about 2-3 Myr, but is older than 
Trapezium. 



4.2. Mass Segregation 

In order to examine the mass segregation in the RCW 36 cluster, we show the stellar 
distribution by splitting the 823 sources detected at Ks with a photometric error of < 0.1 
mag into three magnitude bins (Kg < 12.0, 12.0 < Kg < 14.0, and 14.0 < K5) in Fig. 7. 
The sources of K^ < 12.0 show high central condensation (Fig. 7a). The sources of 12.0 < 
K5 < 14.0 also show central condensation but are less condensed than the sources of K5 < 
12.0 (Fig. 7b). The sources of 14.0 < K^ show no central condensation (Fig. 7c); most of 
these sources are likely to be background stars because their distribution is uniform. These 
figures clearly demonstrate an excess of the brighter stars in the cluster center. This is 
more evident in Fig. 8, which shows the SSND radial profiles for three magnitude bins (the 
same as above). The central position of each radial profiles are the same as Fig. 3. King 
profile fitting for each magnitude bins yields the core radii of Tc = 0.02 pc (Fig. 8a), 0.13 pc 
(Fig. 8b), and 0.24 pc (Fig. 8c, was not fitted well because the surface density decreases in 
the central annulus). 

In young embedded clusters, because most members are considered to be in the pre-main 
sequence (PMS) phase, this tendency could have two possibilities. One is that the sources 
located in the cluster center are young, and the other is that they are massive. The NIR 
excess of the PMS stars is caused by optically thick circumstellar disks/envelopes. Since 
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the circumstellar disks/envelopes become optically thinner with age, the PMS stars which 
possess NIR excess are considered to be younger than those with no NIR excess. Therefore, if 
younger stars concentrate in the cluster center, the distribution of the sources which possess 
NIR excess should be similar to those of the brighter sources. However, as shown in Fig. 8d, 
the radial profile for the sources which possess NIR excess does not show significant central 
condensation. Their King profile fitting yields the core radius of Tc = 0.20 pc which is clearly 
larger than that of the sources of K^ < 12.0. Therefore, the brightness of the sources in the 
central region of the RCW 36 cluster is not due to their youth, but due to their massiveness. 

Mass segregation can be caused 1) by dynamical evolution and equipartition of energy 
after star formation and 2) by the process of star formation itself. If mass segregation is 
caused by dynamical evolution, it occurs in approximately the relaxation time, r^eiaxi of the 
cluster (Bonnell & Davies 1998). The relaxation time is roughly 

O.IAT 

'relax ^^ ~, TT' cross 

InJy 

where N is the number of stars contained in the cluster and Tcross is the crossing time (Binney 
& Tremaine 1987). The crossing time, Tcross, is the time it takes a star to cross the cluster 
{Tcross ~ '^R/vdisp, where R is the radius and Vdisp is the velocity dispersion of the cluster). 
In the case of most open clusters, mass segregation is likely caused by dynamical evolution 
[e.g., Pleiades (Pinfield, Jameson, & Hodgkin 1998)]. However, mass segregation in rich 
and young embedded clusters, such as the Trapezium cluster, is not considered to be caused 
by dynamical evolution because they are too young to have dynamically evolved. Here, we 
estimate the relaxation time for several embedded clusters in table 2. Since the Trapezium 
cluster has a small age and a large number of stars, the relaxation time is significantly larger 
than its age in spite of the small crossing time. In other embedded clusters whose SSND is 
lower, dynamical evolution is generally very unlikely because the cluster sizes are larger, the 
ages are smaller, or both. However, the relaxation time of the RCW 36 cluster is estimated to 
be ~ 2.6 X 10^ yr; this is comparable to the cluster age estimated above. Thus, at least from 
our data, there is a possibility that the mass segregation in the RCW 36 cluster was caused 
by dynamical evolution, because the RCW 36 cluster has smaller radius, moderate number 
of stars, and moderate age. However, this dose not rule out the possibility of primordial 
mass segregation in the RCW 36 cluster because the clusters which dynamical equilibrium 
are established lost all traces of their initial conditions. 
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5. CONCLUSION 

We have presented deep J, iJ, Ks observations of the RCW 36 cluster which hes in a 
C^^O clump in the cloud C of the Vela Molecular Ridge. The main conclusions are as follows. 

1: The RCW 36 cluster has at least ~ 350 members and a cluster radius of ~ 0.5 pc, and 
the radial profile of the stellar surface number density can be well fitted with King profiles 
better than with 1/r profiles. In the central 0.1 pc x 0.1 pc region of the RCW 36 cluster, 
the stellar surface number density is approximately 3000 stars pc~^ This is much higher than 
most of the young clusters within 1 kpc of the Sun. 

2: We measured the NIR excess fraction of the RCW 36 cluster to be 22 %, compared 
with those of other embedded clusters, and estimated the age of the RCW 36 cluster is 
approximately 2-3 Myr. The comparison of background corrected KLF of the RCW 36 
cluster with Trapezium and IC 348 also supports this estimation. 

3: The RCW 36 cluster shows an excess of the brighter stars in its central region. From 
the comparison of the radial profiles with the sources of Ks < 12.0 mag and the sources with 
NIR excess, we conclude that this tendency means that more massive stars are located near 
the cluster center. 

It is pleasure to thank the staff of the South African Astronomical Observatory for 
their kind support during the observations. We also would like to express our thanks to 
the anonymous referee for comments and suggestions, which greatly improved the scientific 
content of the paper. This work is partly supported by Grant-in-Aid for Scientific Research 
of the Ministry of Education, Culture, Sports, Science, and Technology of Japan. 
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Fig. 1. — JHKs three-color composite image of the RCW 36 cluster. (J : hlue^ H : 
green, Ks '■ red) obtained with IRSF/SIRIUS. North is up, and east is to the left. The 
plus represents the peak position of radio continuum (Walsh et al. 1998), and the asterisk 
represents the position of H2O maser (Braz & Epchtein 1982). [Figure separately included 
as gif (to make smaller).] 
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Fig. 2. — Contour map of the stellar surface number density obtained by counting stars in 
a 30" X 30" (~ 0.1 pc x 0.1 pc at 700 pc) area every 15" over the Ks image. The lowest 
contour is 400 stars pc~^, and steps are 200 stars pc~^. Dotted square is the region used to 
determine the completeness limits (see §2). 
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Fig. 3. — Radial profile of the stellar surface number density for the RCW 36 cluster. Stellar 
surface density (in stars per pc^) is measured in circular annuli, centering on the location 
of the peak position of the stellar surface density (see Fig. 2). Profiles are calculated using 
annuli with equal radial steps (=15"). The errors are statistical a/ZV error. The horizonal 
dot-dashed line corresponds to background surface density [level = 2.0 stars in a 30" x 
30" unit cell (see §3.2.3)]. The solid curve represents a King profile fit, and the dashed curve 
represents a 1/r profile fit. 
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Fig. 4. — The color-color diagram for the sources in the RCW 36 cluster with J, H, Ks 
band photometric errors < 0.1 mag. The solid lines represent the loci of dwarfs and giants 
(Tokunaga 1996), the dot-dashed line represents the locus of unreddened Classical T Tauri 
Stars (CTTS) (Meyer, Calvet, & Hillenbrand 1997). The two leftmost parallel dashed lines 
define the reddening band for dwarfs and giants. They are parallel to the reddening vector. 
The tick marks on the reddening line indicate 5 mag intervals in Av from the right edge of 
the CTTS locus. The reddening low of Koornneef (1983), having the slope of 1.7, is adopted. 
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Fig. 5. — The KLFs of the RCW 36 cluster. The bin width is 0.5 mag, and the errors are 
statistical yN error. The solid vertical line represents 90 % completeness limit, a): The 
filled circles with solid line are "raw KLF" (from all sources detected in the Ks band with a 
photometric error of < 0.1 mag), and the open squares with dashed line are the "off-cluster 
KLF" [from the sources detected in the Ks band with a photometric error of < 0.1 mag and 
are located in the off-cluster region (the region whose SSNDs are < background -|- 1 a)]. 
b): The KLF for the RCW 36 cluster members. In order to obtain this KLF, the off-cluster 
KLF (scaled to account for the different areas) is subtracted from the raw KLF. 
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Fig. 6. — The comparisons of the background subtracted KLFs. All magnitudes are shifted 
to absolute magnitudes. The Trapezium and RCW 36 KLFs have been scaled to contain 
the same number of stars as IC 348. The KLFs for IC348 and Trapezium are obtained 
by Muench et al. (2003). The solid vertical line represents our 90 % completeness limit. 
Reddening correction was performed by shifting each KLF toward the brighter magnitude 
by each mean extinction (see §4.1). 
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Fig. 7. — a) - c): The distributions of the sources detected in the K5 band with a photometric 
error of < 0.1 mag [a): the sources of K5 < 12.0, b): the sources of 12.0 < K5 < 14.0, c): 
the sources of 14.0 < K5]. d): The distributions of the sources identified in three bands and 
falling into a region of infrared excess in Fig. 4. 
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Fig. 8. — Radial profiles of the stellar surface number density for the RCW 36 cluster. The 
errors are statistical vN error, a): the sources of Kg < 12.0, b): the sources of 12.0 < 
Ks < 14.0, c): the sources of 14.0 < Kg, d) the sources with NIR excess (see §3.3). The 
photometric errors for all sources are < 0.1 mag. The solid lines represent King profile fits. 
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Table 1. Comparison of the most massive star and the central SSND 



Region 


Distance 


Most massive 


Central SSND 


K 


References 




pc 


Sp type 


stars 


pc ^ 


(Umit) 




RCrA 


150 


B9 


260^ 




16.5 


1 


IC348 


320 


B5 


955^ 




15.0 


2 


NGC2024 


415 


08 


1600^ 




14.0 


3,4,5 


MonR2 


830 


BO 


1900^= 




14.5 


6 


Trapezium 


470 


07 


5600'^ 




17.0 


7 


RCW36 


700 


08-09 


3000^ 




16.3 


this work 



^innner central 0.1 pc 
''yielded by King profile fit 
'^yielded by Gaussian profile fit 

References. — (1) Wilking et al. (1997); (2) Lada & Lada (1995); (3) Lada et al. (1991); (4) 
Bik et al. (2003); (5) Haisch, Lada, & Lada (2000); (6)Carpenter et al. (1997); (7) Hillenbrand 
& Hartmann (1998) 
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Table 2. Estimation of the Relaxation time 



Region 


Size 


^^ star 


"^disp 


^cross 


T-relax 


Cluster Age 


References 




pc 




km s^^ 


[Myr] 


[Myr] 


Myr' 




IC348 


1.0 


300 


1.0^ 


2.0 


10 


2-3 


1,2,3 


NGC2024 


0.88 


309 


1.3^ 


1.3 


7.1 


1 


1,4,5 


MonR2 


1.85 


371 


l.li' 


3.3 


21 


<3 


1,4,6 


Trapezium 


0.24 


780 


2.5= 


0.13 


2.2 


<1 


1,7,8 


RCW36 


0.5 


349 


2.2^ 


0.44 


2.6 


2-3 


this work, 9 



^Av of Ci^O 
^Av of i^CO 



^velocity dispersion of the stars 



References. — (1) Lada & Lada (2003) and references therein; (2) Bachiller, Guilloteau, 
& Kahane (1987); (3) Haisch, Lada, & Lada (2001a); (4) Miesch & Bally (1994); (5) Haisch, 
Lada, & Lada (2000); (6)Carpenter et al. (1997); (7) Jones & Walker (1988); (8)Hillenbrand 
& Hartmann (1998); (9) Yamaguchi et al. (1999) 
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